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We study nuclear effects in the F^{x) structure function in the deep inelastic neutrino reactions on 
iron by using a relativistic framework to describe the nucleon spectral functions in the nucleus. The 

results for the ratio R{x, Q^) = ."'^jv?''^2^ the Gross-Llewellyn Smith(GLS) integral G{x, Q^) = 

dxF-^{x, Q^) in nuclei are discussed and compared with the recent results available in literature 
from theoretical and phenomenological analyses of experimental data. 
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I. INTRODUCTION 

The recent experimental results reported by the NuTeV collaboration [l| on weak charged and neutral current 
induced (anti)neutrino processes on an iron target in the deep inelastic region have emphasized the importance of 
nuclear medium effects. There are many theoretical analyses of the deep inelastic scattering of charged leptons from 
nuclear targets where various nuclear effects like shadowing, anti-shadowing, Fermi motion and binding of the nucleons 
in various kinematic regions have been studied. These have been discussed in several review articles 0, S 13 ; but 
in the case of deep inelastic scattering of (anti)neutrinos from nuclear targets, there are few calculations where the 
dynamical origin of the nuclear medium effects has been studied [1, [E 0, H 0- In some theoretical analyses, nuclear 
medium effects have been phenomenologically described in terms of a few parameters which are determined from 
fitting the experimental data of charged leptons and (anti)neutrino deep inelastic scattering from various nuclear 

■ targets 

\ The differential scattering cross section for the deep inelastic scattering of (anti)ncutrinos from unpolarized nucleons 
in the limit of Icpton mass m; — + 0, is described in terms of three structure functions, Fi{x,Q^), F|'(a;,Q^) and 
F^(x,Q^), where x = ^^-^ — — 2a7i7 Bjorken variable, v and q being the energy and momcntmn transfer of 

leptons. In the asymptotic region of Bjorken scaling i.e. —^ oo, v —> oo, x finite, all the structure functions 
depend only on the Bjorken variable x. In this scaling limit, Fi{x) and i^2^(a;) are related by the Callan-Gross 
relation [l^| leading to only two independent structure functions -Pf (2^) ^^d F^{x) which are determined from the 

I experimental data on deep inelastic scattering of (anti)ncutrinos in the asymptotic region. Scaling violation effects 
have been studied in deep inelastic scattering of (anti)neutrino from nucleon targets using methods of perturbative 
and non-perturbative QCD. The dependence of the structure functions is determined from the evolution given 
by Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations [l^ obtained in perturbative QCD calculations. 
Non perturbative corrections to the structure functions have been also studied by many authors [13, [3 ■ 

Since most of the experimental data on deep inelastic scattering of high energy (anti)ncutrinos are obtained on 
nuclear targets over a wide range of Q^, it is important to study the nuclear medium modification effects on these 
structure functions and F^{x,Q'^) specially in iron where most recent, high statistics, high precision data 

are available The nuclear medium modification effects on the average structure function F^ {x ,Q^)=^{F^^ {x ,Q'^) 
+ F^^{x,Q^)) on isoscalar nuclear targets arise mainly due to modification of the valence quark parton distributions 
in the nuclear medium. The contributions of sea quarks and gluons in a quark parton picture vanish in the case of 
a symmetric sea and arise only when the various sea quark parton distribution are taken to be different. Of much 
interest is the estimate of the size of the nuclear medium effects on the GLS sum rule [19] which is measured over a 
wide range of for deep inelastic scattering of (anti)neutrinos on nuclear targets and has been used to determine 
the QCD coupling constant as [20|. 

In this paper, we study some nuclear medium effects on the nucleon structure function F^{x,Q^) in iron. We use a 
theoretical spectral function to describe the momentum distribution of nucleons in the nucleus. The spectral function 
has been calculated using the Lehmann's representation for the relativistic nucleon propagator and nuclear many 
body theory is used to calculate it for an interacting Fermi sea in nuclear matter. A local density approximation 
is then applied to translate these results to finite nuclei. The method has been earlier used successfully to describe 
the EMC and other nuclear effects in deep inelastic scattering of charged leptons from nuclei [2l[. Our approach 
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describes an alternative method to calculate nuclear medium effects in the nucleon structure functions using nuclear 
many body theory different from the approach used by Kulagin Q. Our model has certain limitations as it only 
attempts to consider the modifications of nucleonic contributions to {x,Q^) arising due to binding energy, off mass 
shell and Fermi motion of the nucleon in the nuclear medium which dominate in the region of x >0.3. In the region of 
0.3> X >0.1, corresponding to the anti-shadowing region, the nuclear medium modification effects on F^{x,Q'^) are 
expected to be small due to vanishing of the pion contribution in this model unlike F^\x,Q'^) where contributions 
of pion and rho mesons play a dominant role [2l|]. We do not consider the shadowing region of 0.0< x <0.1 in this 
paper and that will be considered in the extension of this work. Therefore, the results presented here should be able 
to describe the dominant contribution of nuclear medium effects to F^{x,Q^) in the range of 0.1< x <1. 

In sections 2 and 3, we describe briefly the formalism used to calculate nuclear effects and present our results and 
discussion in section 4. 



II. DEEP INELASTIC NEUTRINO NUCLEON SCATTERING 



The cross section for the reaction 

Mi>i) + N ^l-il+)+X, (1) 



is given by: 



1 2m^ 2M r d^k' 2mi 
^ ^2£;,(k) 2E{v) J 2^3 2£;;(k') ^ > 



i=l lef ^ I jeb \ 1 / 



where / stands for fermions and h for bosons in the flnal state X. The index i is split into I and j for fermions and 
bosons respectively, 

T is the invariant matrix element for the above reaction and is, written as 

^ zT = {^-^ Sz(fc')7"(l - iMk) {qT^^ {X\MN) . (3) 

After performing the phase space integration in Eq.([2]), the double differential scattering cross section evaluated for 
a nucleon target in its rest frame is expressed as. 



<P(^y,u ^ Gf'^ I'k'l / in^ \ T<^P 
dn'dE' (2^)2 |k| W-m^ I '^-'^'''^^ 



L^W^,, (4) 



where Gp is the Fermi coupling constant, mw is the mass of the W boson, l{— e, fi) is a lepton, k is the incoming 
neutrino four momentum and k' is the outgoing four momentum of the lepton, q — k ~ k' is the four momentum 
transfer and U,' ,E' refer to the outgoing lepton. TV is a nucleon, X is a jet of n hadrons consisting of fermions(f) and 
bosons(b) in the flnal state labeled by I and j in the following. 

The lepton tensor for antineutrino(neutrino) scattering L"^ is given by 

^ j^ayp _^ _ J^ ^g^P ± ie'^PP" kpk'„ , (5) 

and the hadronic tensor W^p is defined as 

^5.= is E En / (Sn {}§) n ii) {x\j.\N){x\MNYi2.)HHp+,~±v'd. 

SN X Si i=l ^ ' lef ^ ' ^ jcb \ j J i=l 

where q is the momentum of the virtual VF, sjq the spin of the nucleon and Si the spin of the fermions in X . In the 
case of antineutrino {X\Ja\N) is replaced by {X\,P^\N). 
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The most general form of the hadronic tensor is expressed as |22| : 

where M is the nucleon mass and are the structure functions, which depend on the scalars and p.q. The terms 
depending on Wa^, W5 and Wg in Eq. ([7]) do not contribute to the cross section in Eq. (g]) in the hmit of lepton mass 
mi 0. 

In terms of the Bjorken variables x and y defined as 

r)2 2 P-1 /o^ 

2M V hu M 

we can write the expression for the differential scattering cross section (in the limit of lepton mass m; — > 0) as 



dxdy TT 



xy^F';^'\x,Q^) +[l-y- ^) F::^'\x,Q^) ± xy{l - y I2)F!;^'\x,Q^)^ , 



(9) 



where the + (— ) sign stands for the neutrino (antineutrino) cross section, and the F^''^ {x,Q'^) are dimensionless 
structure functions defined as 

F^^'\x,Q^) - MW^'^'\i,,Q') (10) 
F,^^'\x,Q') = :^W^^'\^,Q') 

In the Bjorken limit of scaling valid in the asymptotic region i.e. — > 00, v —> 00, x finite, the structure functions 
F^'"'^\x,Q'^) are independent of and depend only on the single dimensionless variable x, and satisfy the Callan- 
Gross relation [l^ given as 2xFi{x) = F2{x). Using this relation, the cross section in Eq.® is described in terms of 
two independent structure functions -^2(2^) and F^^x). In the quark parton model of deep inelastic scattering, in the 
Bjorken scaling limit, these structure functions are determined in terms of parton distribution functions for quarks 
and antiquarks. 

Specifically, the structure function F3{x) is given as: 



F^P = 2[d{x) + s{x) - u{x) - c{x)] , 
Fp = 2[u{x) + s{x) - d{x) ~ c{x)] , 
Fl^^ = 2[u{x) + c{x) - d{x) - s{x)] , 



Fg^" = 2[d{x) + c{x) - u{x) - s{x)] , 
The average structure function F^ [x) on isoscalar nucleon target defined as 

F^{x)^\{F^,^ + Fr) 

is given by 

F^^ [x) = [uy{x) + dv{x) + s{x) — s{x) + c{x) — c{x)], 

where Uv{x) — u{x) — u{x) and dv{x) = d{x) — d{x) are the valence quark parton distributions. 

Thus, for an isoscalar target and a symmetric sea, F-^ (x) structure function is given in terms of valence quarks u„ 
and dy which satisfies the Gross-Llewellyn Smith sum rule : 

/ Fi^{x)dx = 3. (11) 
^0 
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FIG. 1: Self-energy diagram of the neutrino in the nuclear medium associated with the process of deep inelastic neutrino-nucleon 
scattering. The imaginary part is calculated by cutting along the horizontal line and applying the Cutkosky rule for putting 
the particles on the mass shell. 



In the non- asymptotic region, scaling violations occur and the structure function Fi (x) acquire dependence which 
are calculated with the DGLAP equations of evolution obtained using perturbative QCD. As a consequence, the 
Callan Gross relation and the Gross-Llewellyn Smith sum rule are modified as follows: 

2.F,(x,Q^)=F.(.,Q^)^^^^, (12) 
where R=f^ is the ratio of the cross section of longitudinally to transversely polarized W bosons [23*] and 

S{Q') = F,M')dx = i{l-^-a{nf)[^)" -b{n,)[^)'Y^^ (13) 

where a{nf) and bjuf) are functions of the number of quark flavors accessible at a given and Aht is the higher 
twist correction 1241. 



III. NUCLEAR EFFECTS IN NEUTRINO SCATTERING 



When the reaction given by Eq.(IT]) takes place on a nucleon in the nucleus, nuclear effects have to be considered. 
There are two main nuclear effects. Firstly, a kinematic effect which arises as the struck nucleon is not at rest but 
is moving with a Fermi momentum in the rest frame of the nucleus, leading to a Lorentz contraction of the incident 
flux used in deriving Eq. ([¥]). Secondly, the more important dynamic effects which arise due to Fermi motion, Pauli 
blocking and strong interaction of the initial nucleon in the nuclear medium. 

In a nuclear medium the expression for the cross section given in Eq.Q is modifled as: 

dV,^, _ Gl |k'| / m'^ \\c.0,,,A 



d^l'dE' (2^)2 |k| W-m^^J ""-'^"'"/S' 

where W^p is the nuclear hadronic tensor defined in terms of nuclear hadronic structure functions WiAix, Q^) through 
Eq.©. Accordingly, the dimensionless structure functions Ff^[x,Q^) are defined in terms of Wf-{x,Q'^) through 
Eq. ()10ll . In our formalism the neutrino nuclear cross sections are obtained in terms of neutrino self energy S(fc) in 
the nuclear medium which also defines the dimensionless nuclear structure functions F^(x,Q^). A comparison with 
Eg. ([36)1 gives the nuclear hadronic tensor W^^. With proper choice of tensor components a,/?, the dimensionless 

nuclear structure functions Ff'{x,Q'^) are obtained [2l|. 

The neutrino self-energy in nuclear matter corresponding to Fig.l is given by. 



Gf 4: f d^k' 1 ( raw 



m= (-*)^— 7^b4 772 ^ L^pli'^^iq), (15) 
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where Lafj is given by Eq.([5]) and U"^{q) is the W self-energy in the nuclear medium and is written with the help of 
Fig. 2 as: 



x{X\J'^\N){X\J^\Nn2nr6\q + p-j:Up'^. (16) 

In the above expression Gi{p[) and Dj{p'j) are respectively the nucleon and meson relativistic propagators in the final 
state which are taken as the standard free relativistic propagators [25|. G{p) is the nucleon propagator with mass M 
and energy E(jp) in the initial state, which is calculated for a relativistic nucleon in the interacting Fermi sea. 

For this we start with the relativistic Dirac propagator G(p) for a free nucleon, which is written in terms of the 
contribution from the positive and negative energy components of the nucleon described by the Dirac spinors u(p) 
and w(p) using their appropriate normalisations '25'| as 

p2 _ M2 E{p) \ pO - E{p) +i€ + £;(p) -ie j ^ ^ 

The nucleon propagator G(p) is then calculated by making a perturbative expansion of G{p) in terms of G'^(p) 
given in Eq. ()17p by retaining the positive energy contributions only (the negative energy components are suppressed). 
This perturbative expansion is summed in ladder approximation to give [2l|: 

T\/T 1 

M ^ Ur{v)ur{v) sr^ M ^ Us{Vi)Us{p) 

+ ij(p)Z. pO_^(p) l^^P 'P^i?(p) P^-E{p) ^ 

- Ur{v)Ur{-p) 

- i?(p) 4- pO _ ^(p) _ ^,,(p) 2A'(pO^ p)^.,.(p)4§y 

where Yi^ (p^ ,p) is the nucleon self energy in nuclear matter taken from Ref. [26| . The relativistic nucleon propagator 
G(p) in a nuclear medium is then expressed as 21[: 



-— > M^pu^p / du— -+ / duj^ — — 



(19) 



where Sh{uj,p) and Sp{uj,p) being the hole and particle spectral functions respectively, which are given by [2ll [26|: 
^'^'^'P^ = ^(p" - £;(p) - ^i?eS^(p0^p))2 + (_M.j^,s^'(pO,p))2 

for p^ < ]i 

p. 1 ^^-^"(^°'^) .21) 

' ^(pO-i?(p)-4|)i?eS^(pO,p))2 + (^/mS^(pO,p))2 ^ ^ 

for > /z. 

The normalisation of this spectral function is obtained by imposing the baryon number conservation following 
the method of Frankfurt and Strikman [27]. For this purpose we evaluate the electromagnetic form factor at q=0, 
assuming baryons have unit charge for the purpose of normalisation, corresponding to Fig. 2a, i.e. 

{N\B>'\N) = u{p)ru{p) - 5 = 1, p'' = {E{p),p) (22) 

When the nucleons are in the nuclear medium the many-body diagram shown by Fig. 2b is evaluated [2l| and gives 

{A\B^\A) = ~j ^yzrr[G(/,p)7le'^^. (23) 
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FIG. 2: Electromagnetic form factors for the case (a) free nucleon, (b) Fermi sea with B baryons. 



where V is the volume of the normalisation box and exp{ip'^rj), with rj 0^, is the convergence factor for loops 
appearing at equal times. 

Using the expression given by Eq. p9p , it can be seen that the convergence factor limits the contribution of the hole 
spectral function and this gives 



onshell 



2V 



(27r)3 E{p) 
(27r)3 £:(p) M 



Sh{uj,p)duj 
Sh{uj,p)duj 

-OO 

S,,iLO,p)dLU = B-f. 

-oo 



2M 



(24) 



It is to be noted that in the last step we have imposed that this matrix element gives the right current with B baryons, 
in analogy to the expression given by Eq. (P^ . is the four momentum of the nucleus. The operator {/jj + M) 
comes from ^Ur(p)ur(p) which depends only on p, and that corresponds to a free particle with {p^ — E{p),p), 
therefore, the operator (^+ M) is written on shell. The last step in Eq. ([M)) ) is obtained after performing the trace. 
Evaluating Eq. in the rest frame of the nucleus we obtain the normalisation as 



2V 



d^p 
(2^ 



Sh{uj,p)duj = B 



(25) 



In the local density approximation, we do not have a box of constant density, and the reaction takes place at a point 
r, lying inside a volume element d^r with local density Pp(r) and Pn(r) corresponding to the proton and neutron. 
Therefore, the upper limit in the integration over nucleon momentum in Eq. (|25p is the local Fermi momentum kp ^ (r) 
of the nucleon given by: 



1/3 



(26) 



This makes the spectral function Sh{uj,p) density dependent i.e. Sh{u!,p, kpij)) and the normalisation condition given 
in Eq. (P5)l is modified to 



_fp_ 

(2^ 



Shiuj.p, kFj,^„{r))duj = Pp,nir) 



(27) 



For a symmetric nuclear matter of density p(r), there is a unique Fermi momentum given by kp{r) = [|7r^p(r)] 
for which we obtain 



1/3 



d^p 
(2^ 



Sh{iL>,p, kF{r))duj = p(r), 



(28) 
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where p{r) is the baryon density for the nucleus. 

This leads to the normalisation condition [2^ given by 

d^^j^f Sh{oJ,p.p{r))<ko = A (29) 

In the antineutrino case the expressions obtained are very similar. L'^^ appears, as in Eq. ([5]), with a minus sign in 
front and in the W self-energy, Eq. (fT6|) . we have {X\Jl\N), instead of {X\Ja\N). 

The probability per unit time for the neutrino to collide with nucleons when traveling through nuclear matter is: 

r(A:) = -;|^ImE(fc), (30) 

and the cross section for neutrino scattering from an element of volume (fir and surface dS in the nucleus is given by 

da = TdtdS = r4 dldS = -d\^T d\ = Im E d\ . (31) 

dl V \k\ \k\ 

Using Eq. ([3T|) in Eq. (fT5|) , we get the expression for the total scattering cross section in the local density approximation 
as 



The imaginary part of the neutrino self energy in Eq. (|32p is evaluated by means of the Cutkosky rules [25[ by 
cutting the Feynman diagram shown in Fig. 3 along the dotted line which puts the particles corresponding to the cut 
propagators on the mass shell by replacing the fermion and meson propagators by their imaginary parts as 

Y,{k) 2iImT.{k) 
D{p/) ^ 2i9{po,) ImD{p/) 
G{pi') ^ 2ie{poi') ImG{pi') 

> 2TrS{k'^ -mf). (33) 
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After performing the fco', po' and po integrations for all momenta in Eq. (|5^ using Eqs. p^ and we get the 

differential scattering cross section which is written in the local density approximation as: 



dn'dE' (2^)2 |k| W-ml,J 

where 



W^f,^A j d\ j ^ J'jpoJL.s^^p0^p^p^^))w^^^,q), . (35) 

with W^^(p, q) is given by Eq.®. 

Note that the factor in Eq. ([55]) comes naturally in our formalism, when we perform the various momentum 
integrations in Eq. (|32p to calculate the imaginary part using Eqs. (|19p and ([55]) . This can be physically understood 
as a kinematic factor which appears in the cross section defined in Eq. ([2]) for a nucleon moving with momentum 
p = {E, p) in the rest frame of the nucleus leading to the equation: 

d^'dE' (27r)2 |k| ^(p) U'-"'4/ ''''' "/^^ ^ ' 

When this cross section for a nucleon target moving with momentum p^ in the rest frame of the nucleus of density 
pir) and weighted with the spectral function Ship*^ ,p, p{r))., is summed over all the nucleons in the nucleus, it leads 
to Eq. p4|) . Similar equations are also obtained in Refs. [3.. ,9.. ,29.. ,3d| . 
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FIG. 3: Self-energy diagram of the W-boson in the nuclear medium. The imaginary part is calculated by cutting along the 
horizontal line and applying the Cutkosky rules when putting the particle on mass shell. 

Now to evaluate F^, we calculate the components xy on both sides of Eq. (pSj). We get by taking q along the z 
axis: 



2Ma 



and for the right hand side we have for the nucleon moving with momentum p: 



M2 



2M^ 



(37) 



(38) 



Since we have 



(39) 



(40) 



with X as defined in Eq. ([8]) and xat is the Bjorken variable expressed in terms of the nucleon variables , {p^ , p), in 
the nucleus 



Xn 



2p.q 2{pQqo - pzq^)' 



we obtain the expression for Fj^{x) in the Bjorken limit 



d^p M r . Oo / / ^^ 
dp^ShiP^.p^pir)) 



xn 



Mq, 



Fi'{xN,Q'), 



where the contribution of W2 appearing in Eq. (j38p vanishes after momentum integration. Defining 7 as 



we get 



F^{x,Q^) = 4 / 



qz 

'^7 



d^p 
(27r)3 



4A^2.x2^^/' 



dp'Sf,ip'\p, p{r))F{p, Q')Fi' {XN,Q') 



(41) 



(42) 



(43) 



(44) 



where F(p, Q^) = -J^ , -, ^ 

This is our main equation which describes the modification on F^[x,Q^) due to nuclear medium effects and was 
earlier obtained in Rcf. Idl. 
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FIG. 5: Results for the ratio R = at = 5 GeV^ by different autliors. Solid line: this work using MRST2004 

NNLO parton distribution function; double dashed-dotted line: Hirai et al. short dashed line: NuTeV collaboration [ij, 

dashed-double dotted line the results of Kulagin and the results of Kulagin and Petti [p| shown by dashed-dotted line. 



IV. RESULTS AND DISCUSSION 



We obtain the numerical results for F^{x^ Q^) using Eq. p4|) and Su given by Eq. ((20|). The density for ^^Fe nucleus 
is taken to be a two Fermi parameter distribution i.e. 

p(r) = (45) 
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■ Barabash et al.[36] 




Q^eV^) 

FIG. 6: AGLS^^{3~ Jg F3^{x,Q^)dx) vs 



where the parameters cl=4.106fni and c2=0.519fm and are taken from Ref [3l|. While using the spectral function 
Sh{uj,p, p{r)) in Eq.([33]), we ensure that it gives the correct binding energy for the iron nucleus by varying the free 
piece of the real part of the nucleon self energy in Eq. ([^01) . 

This is done by calculating the average kinetic and total nucleon energy given by: 



< T >= - Id 



(46) 



<E>=^jd\ Su{p\p,p{r)) p'dp' , (47) 

and the binding energy per nucleon given by [2l|: 

\Ea\ = -1{<E-M>+^^<T>) (48) 

After the free parameter in is fixed, the expected values of the kinetic energy < T > and the total nucleon energy 
< E > are obtained as 30MeV and 48MeV for the case of ^^Fe which are in good agreement with other models and 
experiment (See e.g. Table II of Ref. [3^ ) . A similar agreement is also obtained with other nuclei. The structure 
function F-^( x, Q^) is numerically evaluated. For F^{xn,Q'^) we use the MRST2004 NNLO parton distribution 
functions [33]. The evolution of F^{x,Q^) is assumed to be given by the evolution of GLS integral S{Q^) 
given in Eq. (12). The strong coupling constant as{Q'^) is calculated using the variable flavour number evolution 
equation with A — 251MeV for nj — A and A — l78MeV for Uf — 5 [24l| which gives 0:5(^2)=. 1153. The constants 
a{nf) and b{nf) in Eq. (12) are taken from Larin and Vermarseren |24| . The target mass corrections are incorporated 
using the prescription of Petti and Kulagin Q following the work of Georgi and Politzer [33| . 

In Fig. 4, we show the results for the ratio K{x,Q'^)— ^^J^^^^^l^ as a function of x for various values of Q^. For 

Q^ = 10GeV^ , we get a suppression for x < a;mm=0.7, beyond which we get an enhancement. This value of Xmin 
decreases with increasing Q^. For x > x,nin, the ratio increases very fast and becomes larger than unity a.s x —> 1. 
This is mainly due to the Fermi motion of the nucleons. In the region of 0.3< x <1, this behavior is very similar to 
that seen in the EMC effect with charged lepton and (anti)neutrino scattering from the nuclear targets in the analyses 

of the structure function F^'^x, Q^) and F^''^'"'^\x, Q^)- 
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FIG. 7: The contribution to the GLS integral from the region x=0.5 to x=l for different using MRST NNLO[33| parton 
distribution function(PDF) and its comparison from the results of Kim et al.[20(]. The dashed line is the result with the free 
case, solid line is the result with nuclear medium effects obtained in our model with spectral function including QCD and TMC 
corrections, and the solid square with error bars are the results of Kim et al.[20| 



1.1 




FIG. 8: The x dependence of the nucleus/nucleon ratio of the GLS integral from the nuclear effects to the free case i.e. the 
'■'^^^^ ji"'" fW(^)^^ at = 5Ge\/^ using MRST NNLO[33| parton distribution function(PDF). The dashed line is the result 
of Kulagin . 



In the region of 0.1 < x <0.3, we do not see any enhancement due to the anti-shadowing effect in F^(x^ Q^) as found 
in F^''(a;, Q^) [2l[. In the present calculation, this is due to the neglect of the mesonic contributions, responsible for 
the enhancement in F^'^ix, Q^). However, for F^(x, Q^) these mesonic contributions are expected to be small because 

of the vanishing of the dominant term due to pions. Even in the case of F2'^^'^\x, Q^), the recent [Tsj estimates of 
nuclear medium effects in the anti-shadowing region, do not see any enhancement in this region, in disagreement with 
the other results available from earlier analyses [lOi [O, ■ This is not unexpected as these structure functions, at 
low X, probe different contributions of the parton distribution functions in the nuclear medium. This suggests that 
the nuclear medium modification effects could be different in F^'^x, Q^), F^''^''-'^\x, Q^) and F^''^'"'^\x, Q^), and calls 
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for a comprehensive study of nuclear medium modification effects on these structure functions specially in the region 
of small X. 

In Fig. 5, we compare our results for 'R{x^Q^) with the results of Tzanov et al. [l|, Kulagin and Petti [s'l, Kulag in i 
and Hirai et al. [ifl]. While the work of Kulagin ^5,] and Kulagin and Petti [1, Q use a nuclear model to calculate the 
nuclear effects which shows a dependence, the work of Tzanov et al. [l| and Hirai et al. [l3| are phenomenological 
analyses, which assume the nuclear effects to be independent of . In this figure, the results for R(x,Q^) for 
= 5GeV^, are presented. We find a suppression in F^{x,Q'^) for x <0.7 and an enhancement thereafter, which 
are respectively smaller than the results of Ref. Q, but are larger than the recent results of Kulagin and Petti 0]. 
It should be noted that these latter results 0] give suppression in the region of 0.4< x <0.8 and enhancement for 
X >0.8, which are smaller than the present results and the results obtained earlier in Ref. 01 ■ When compared with 
the results of Tzanov et al. and Hirai et al. T^, we find a smaller suppression in the region 0.5< x <0.7. In the 
region 0.7< x <0.8, we find an enhancement while they obtain a suppression. The results of the phenomenological 
analyses shown in Fig. 5 [l|, [l0| use the experimental data on F^'\x, Q^) in charged lepton scattering to estimate 
the nuclear effects. It can, therefore, be inferred that we find nuclear effects in the case of F^lx^Q"^) different than 
obtained in the case of F2 '''{x, Q^). This is true for all in our model. 

The GLS sum rule fEq fTTI) provides a benchmark to test various models used for the calculation of F^{x, Q^). We 
find that in the —> 00 limit, the correction to the GLS sum rule in our model comes from the off shell modifications 
to the nucleon spectral function. It is easy to see this from Eq. . The nuclear GLS sum rule in the limit ^ 00 
(i.e. 7 ^1) is given by: 

f^F^{x)dx^A Jd\ J'^ dp^Sn{p\p.p{r)){?^^^F^{x^)dx^ (49) 



In the limit of noninteracting nucleons this trivially reproduces the GLS sum rule for free nucleons. 

The dependent nuclear effects in the GLS integral enter through the factor 7 and the dependence of 
F^ {xjsi, Q^) in Ea. (|151) when it is integrated over x. The dependent nuclear corrections to the GLS sum rule are 
thus linked to the dependent perturbative and non-perturbative QCD effects appearing in F^ {xn, Q^)- In Fig. 6, 

we show the Q-^ dependence of the nuclear effects of the GLS integral, where we plot AGLS=i(3 — F ^{x, Q^)dx) 
12 ' ' 



as a function of Q^. The experimental results from CCFR collaborations [20|, CHARM collaborations j35j and IHEP- 
JINR collaborations [3^ are also shown. The behavior of AGLS has been found to be in agreement with the present 
available experimental results. In this figure, we have also shown the theoretical results obtained by Qiu and Vitev [3]. 
Our results are in agreement with the results of Qiu and Vitev 7] for > 5GeV^ where theoretically the suppression 
is found to be larger than the experimental results. For < hGeV^ , we find a larger suppression compared to the 
central value of the experimental result and both theoretical values are within the experimental errors. The nuclear 
corrections to the GLS sum rule found by Kulagin and Petti are quite small due to the cancellation of nuclear 
shadowing, not included in our model, and off shell effects. 

We show in Fig. 7, the value of S{x,Q^)= dxF^{x, Q^) for x=0.5, with and without the nuclear effects and compare 
them with the experimental results of Kim et al.[20]. The nuclear effects tend to reduce S(x,Q^) but the reductions 
are smaller than the experimental uncertainties. In Fig. 8, we show the results for S(x,Q^)=J dxF^{x,Q^) vs x 
at Q^=5GeV^ , where we have also shown the result obtained by Kulagin Q. We find that nuclear medium effects 
leading to a suppression in 0.2< x <0.55 region and an enhancement in x >0.55 region are smaller than obtained by 
Kulagin 

To summarize our results, we have studied nuclear effects in the structure function Fj^(x,Q^) in the iron nucleus 
using a many body theory to describe the spectral function of the nucleon in the nuclear medium for all Q^. The 
nuclear effects are found to decrease the value of the structure function for x < Xrn,in=0-7 and increase at higher 
X > Xmin- The parameter Xmin is found to be dependent which decreases with Q^. The results are compared 
with other theoretical analysis of Kulagin ,5,] and phenomenological analyses of Tzanov et al. ^Ij and Hirai et al. p^ . 
The effect of the nuclear medium modifications on the Gross-Llewellyn Smith [l^l sum rule and its Q^-dependence 
has also been studied. In general nuclear medium effects decrease the value of GLS integral for all . 
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